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We present critical fields, thermally-activated flux flow (TAFF) and critical current density of 
tetragonal phase /3-FeSe single crystals. The upper critical fields H C 2(T) for H||(101) and H_L(101) 
are nearly isotropic and are likely governed by Pauli limiting process. The obtained large Ginzburg- 
Landau parameter k ~ 72.3(2) indicates that /3-FeSe is a type-II superconductor with smaller pen- 
etration depth than in Fe(Te,Se). The resistivity below T c follows Arrhenius TAFF behavior. For 
both field directions below 30 kOe single vortex pinning is dominant whereas collective creep be- 
comes important above 30 kOe. The critical current density J c from M-H loops for H||(101) is about 
five times larger than for H_L(101), yet much smaller than in other iron-based superconductors. 

PACS numbers: 74.70.Xa, 74.25.Op, 74.25.Wx, 74.25.Sv 
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I. INTRODUCTION 

Among iron-based superconductors, the tetragonal /3- 
FeSe has trigged great interest because of simplest struc- 
ture and superconductivity below about 8 K without any 
carrier doping. 1 Besides /3-FeSe, the superconductivity 
has also been discovered in binary Fe(Te,Se) and Fe(Te,S) 
materials These materials also have similar Fermi sur- 
face to FePn (Pn = P, As)-based superconductors^ even 
though they have FeCh (Ch = Te, Se, S) layers stacked 
along the c-axis as opposed to FePn layers. Understand- 
ing this similarity is rather important. Binary FeCh ma- 
terial, /3-FeSe has another notable characteristics. Appli- 
cation of pressure leads to a significant enhancement of 
T c up to 37 K at around 9 GPa, the third highest known 
critical temperature for any binary compound. - 

In order to study anisotropic and intrinsic physi- 
cal properties of materials, single crystals are required. 
When compared to Fe(Te,Se) and Fe(Te,S), extremely 
complex phase diagram of FeSe makes impurities such 
as a-FeSe and FeySes ubiquitous in as-grown crystals, or 
sometimes polycrystals impeding the understanding of 
/3-FeSe£-i£ 

Here, we report intrinsic superconducting properties 
of /3-FeSe single crystals. These include critical fields 
H C 2 and H&, thermally-activated flux flow (TAFF) be- 
havior, and the critical current density J c . Our results 
show that /3-FeSe is a type-II superconductor with large 
Ginzburg-Landau parameter k and smaller critical cur- 
rent density when compared to other iron-based super- 
conductors. Single vortex pinning dominates vortex dy- 
namics below 30 kOe, whereas collective creep becomes 
important at higher magnetic fields. 



II. EXPERIMENT 

Details of synthesis and structural characterization are 
explained elsewhere^ Thin Pt wires were attached to 
electrical contacts made of Epotek H20E silver epoxy for 



a standard four-probe measurement, with current flow- 
ing in crystal plane. Sample geometry were measured 
with an optical microscope Nikon SMZ-800 with 10 pm 
resolution. Magnetization and resistivity measurements 
were carried out in Quantum Design MPMS and PPMS, 
respectively. 



III. RESULTS AND DISCUSSION 

Figure 1(a) shows the temperature dependence of the 
in-plane resistivity p(T) below 300 K. The residual resis- 
tivity ratio (RRR) of 14 is double that of the hexagonal 
shape crystals, 6 indicating good sample quality. The cur- 
vature of p(T) changes at about 100 K (Fig. 1(a) inset 
(a)) due to structural and magnetic transitions in agree- 
ment with previous results fi£J£ With further decrease 
in temperature, superconductivity emerges with T Ct0nset 
~ 11.4 K and T Cj0 ~ 8.1 K (Fig. 1(a) inset (b)), similarly 
to reported values in the literature^ Fig. 1(b) presents 
the ac susceptibility of FeSe single crystal for field per- 
pendicular to the crystal plane. The corresponding su- 
perconducting volume fraction at T = 1.8 K is about 70 
%, confirming the bulk nature of superconductivity. The 
XRD pattern of a single crystal (Fig. 1(b) inset) reveals 
that the crystal surface is normal to (101) direction, also 
similar to previous results 

Figure 2 (a) and (b) show temperature dependence of 
p{T) in various fields for H||(101) and H_L(101). The 
T c shifts to lower temperature without obvious broad- 
ening for both directions with the increase in magnetic 
fields. Temperature dependence of the upper critical 
fields H c2 (T) was determined from the resistivity drops 
to 90%, 50%, and 10% of the normal-state resistivity 
p n (T,H) (Fig. 2(c)). The normal-state resistivity was de- 
termined by linearly extrapolating the normal-state be- 
havior above the onset of superconductivity. The H C 2 (T) 
curves are nearly linear and the initial slopes dH^/ oIT\t c 
are given in Table 1. The slopes are nearly identical for 
both field directions and the anisotropy of the upper crit- 
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FIG. 1. (a) Temperature dependence of the in-plane resistiv- 
ity p(T) of /3-FeSe single crystals. Inset (a) shows derivative 
of resistivity data dp/dT as a function of temperature. Inset 
(b) shows enlarged resistivity curve near T c . (b) Tempera- 
ture dependence of dc magnetic susceptibility of /3-FeSe single 
crystals for H = 10 Oe along the (101) plane. Inset: single 
crystal XRD pattern of /3-FeSe. 



ical field j{T) = H c2MMW1) (T)/H c2 ,h\\(ioi) i T ) is also 
isotropic within the experimental error. 

Within the weak coupling BCS theory^ and using the 
slope determined from the midpoint of resistive tran- 
sition with T c — 9.8 K we can estimate £f C 2(0) = - 
0.693T c (dff c2 /dT| T J = 180(4) kOe for both field di- 
rections. The results are close to the Pauli paramag- 
netic limit H P (0) = 1.84T C = 180(4) kOe. 14 It implies 
that the spin-paramagnetic effect may be the dominant 
pair-breaking mechanism in FeSe for both field direc- 
tions, similar to Fe(Te,Se) and Fe(Te,S)J£J£ The su- 
perconducting coherence length £(0) estimated using the 
Ginzburg-Landau formula if C 2(0) = $o/27r£ 2 (0), where 
$o = 2.07xl0~ 15 Wb is the flux quantum, is f(0) = 
4.28(5) nm, which is somewhat larger than Fe(Te,Se) and 
Fe(Te,S)^^ 

According to the thermally-activated flux flow (TAFF) 
theory, the Inp — l/T in TAFF region can be described 
using Arrhenius relation ) 17 i 18 



lnp(T,H) = lnp (H) - U (H)/T 



(1) 



where lnpo(H) = Inpo t +Uq(H) /T c is the temperature- 
independent constant and Uq(H) is the apparent acti- 
vated energy. Hence, the lnp(T,H) vs. l/T should be 
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FIG. 2. Temperature dependence of the resistivity p(T) of 
/?-FeSe single crystals for (a) H(101) and (b) H(101) in the 
magnetic field up to 90 kOe. (c) Temperature dependence of 
the resistive upper critical field H C 2{T) corresponding to three 
defined temperatures for both field directions (see text). 



linear in TAFF region. As shown in Fig. 3 (a) and (b), 
the Arrhenius relation (solid lines) can fit the experimen- 
tal data very well for both field directions. The results 
are shown in the common logarithmic scale in the fig- 
ures, but we calculate them in the natural one. The 
obtained Uq are similar for both field directions. They 
are comparable to that in Fe(Te,S) and much smaller 
than in Fe(Te,Se)J^^S The good linear behavior indi- 
cates that the temperature dependence of thermally ac- 
tivated energy (TAE) U(T, H) is approximately linear, 
i.e., U{T,H) = U a (H)(l - T/T C )M^. The log p{T,H) 
lines for different fields extrapolate to the same temper- 
ature T cross , which should equal to T C J£ The extrapo- 
lated temperatures are about 11.1 K for both H|| (101) 
and H_L(101). Moreover, lnpo(H) — Uq(H) show linear 
behavior for both field directions (Fig. 3(c)). Fits using 
lnpo(H) — Inpof + Uq(H)/T c , yielded values of po/ an d 
T c 37(1) mn ■ cm and 11.2(1) K for H||(101) and 39(2) 
mfl ■ cm and 11.2(1) K for H_L(101). The T c values are 
consistent with the values of T cross within the error bars. 
The U (H) shows a power law (U (H) - H~ a ) field de- 
pendence for both directions (Fig. 3(d)). For H||(101), 
a = 0.25(6) for H < 30 kOe and a = 0.68(6) for H > 30 
kOe; For H_L(101), a = 0.26(2) for H < 30 kOe and a = 
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TABLE I. Superconducting parameters of /3-FeSe. 





{dH c2 /dT) Tc 


H c 2 } mid(0) 


C(0) fl>(0) Hci(0) 


A(0) H c (0) k(0) 


(K) 


(kOe/K) 


(kOe) 


(nm) (kOe) (Oe) 


(nm) (kOe) 




Onset Middle Zero 








H||(101) 9.8(2) 


-30.3(6) -26.5(6) -25.4(4) 


180(4) 


4.28(5) 180(4) 75(2) 


309(4) 1.76(3) 72.3(2) 


H-L(lOl) 9.8(2) 


-29.6(6) -26.6(6) -25.5(4) 


180(4) 


4.28(5) 180(4) 





(a) 



§ 
& -2 



-4 
0.08 



^ ° 0T 


0.5T 


\^ ° 1T 


3T 


Is^x 5T 


° 7T 








H//(101) 







(b) 





o 0T 


0.5T 


V 


o 1T 


3T 




o 5T 


7T 












Hl(101) 









0.12 0.16 0.20 0.24 0.28 
1/T(1/K) 



2 

1 


-1 
-2 
-3 

0.08 0.12 0.16 0.20 0.24 0.28 
1/T (1/K) 




FIG. 3. Logp(T,H) vs. 1/T in various field for (a) H||(101) 
and (b) H±(101). The corresponding solid lines are fitting 
results from the Arrhenius relation, (c) lnpo(H) vs. Uo(H) 
derived from Arrhenius relation for both field directions. The 
solid lines are linear fitting results, (d) Field dependence of 
Uo(H). The solid lines are power-law fitting using Uo(H) ~ 
H~ a . 



0.70(9) for H > 30 kOe. The weak power law decreases of 
Uq(H) in low fields for both field directions implies that 
single- vortex pinning dominates in this region^ followed 
by a quicker decrease of Uq(H) in field which could be 
related to a crossover to a collective flux creep regime^ 

Low field M(H) at various temperatures for H||(101) 
and H_L(101) are shown in Fig. 4(a) and (b), respec- 
tively. All curves exhibit linear behavior for low fields 
and then deviate from linearity at different field values 
for different temperatures. The values of H c \ are de- 
termined by examining the point of deviation from the 
linear slope of the magnetization curve. The temperature 
dependence of H c i (T) for both field directions are shown 
in Fig. 4(c). For H||(101), the H cl (T) can be well fitted 
using the formula H cl (T) = H cl (0)[l - {T/T c ) 2 } and the 
obtained H, 



ci,H||(ioi)(0) is 75(2) Oe. For BLL(lOl), it is 
difficult to estimate the H cl H j_^ 101 )(T) and obtain reli- 
able fits. This is due to the small obtained value for H cl 
with relative large error and significant demagnetization 




FIG. 4. Low field M(H) of /3-FeSe single crystals at various 
temperatures for (a) H||(101) and H_L(101). (c) Temperature 
dependence of H c i(T) for both field directions. The solid 
line is the fitted lines using H cl (T) = H cl (0)[l - {T/T c ) 2 ] for 
H||(101). 



factors. 

Since Fe(Te,Se) and Fe(Te,S) superconductors are in 
the dirty limit, we assume the same for FeSej22, Using esti- 
mated values of /f c i ) H-||(ioi)(0), H C 2(0) and £(0), we eval- 
uate additional parameters using expressions H C 2(0) = 
V2kH c (0) and H cl (0) = ^^(lnK+0.08) where k = A/£ 
is the Ginzburg-Landau (GL) parameter and H c (0) is the 
thermodynamic upper critical field at T = K.— We ob- 
tain K miol) (0) = 72.3(2), ff c , H ||(ioi)(0) = 1.76(3) kOe, 
and penetration depth for H|| (101) A ff || (1O i)(0) = 309(4) 
nm which is somewhat smaller than in Fe(Te,Se)i22 All 
superconducting parameters are listed in Table 1. 

Figures 5(a) and (b) show isothermal magnetization 
curves M(H) at various temperatures for H||(101) and 
H_L(101). The shape of the M(H) curves confirms that 
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FIG. 5. M(H) loops of /3-FeSe single crystals at various 
temperatures with field up to 50 kOe for (a) H||(101) and 
H_L(101). (c) In-plane and (d) interplane superconducting 
critical currents as determined from magnetization measure- 
ments using the Bean model. 



/3-FeSe is a typical type-II superconductor. The data 
exhibit a central peak at zero magnetic field and then 
the magnetization decreases continuously with increas- 
ing magnetic field. On the other hand, there is a 
weak ferromagnetic (WFM) background superimposed 
on the superconducting M(H) curve for H||(101). This 
WFM background possibly arises due to vacancy-induced 
magnetic cluster. 11 From the irreversible parts of the 
M(H) loops, the critical current can be determined us- 
ing the Bean modeli^^ For a rectangular-shaped crys- 
tal with dimension c < a < b, when H_L(101), the in- 



plane critical current density J, 



ff-L(101) 



(H) is given by 



Jc HW1) {H) = 20AM(H)/[a(l - a/36)], where a and 
b (a < b) are the in-plane sample sizes, AM(H) is 
the width of the magnetic hysteresis loop in emu/cm 3 . 
On the other hand, for H||(101), there are two different 
current densities: the vortex motion across the planes, 
J^ ross {H) and that parallel to the planes, JP ara {H). As- 



suming a,b > (c/3) • JP ara (H)/J^ ross (H)^. we obtain 
Jf oss (H) fa 20AAf (H)/c. The magnetic field depen- 
dence of J^ ross (H) and Jc UW1) {H) is shown in Fig. 
4(c) and 4(d), respectively. It should be noted that 
for H||(101), the WFM background has only minor ef- 
fect on the calculation of AM(H), because of very weak 
moment contribution of WFM. The error of calculated 
J^ ross {H) is about 10% without subtracting WFM back- 
ground. It can be seen that the J^ ross (0) and J^^ 101 ' (0) 
at 1.8 K are about 2.2xl0 4 and 4xl0 3 A/cm 2 , which are 
much smaller than in Fe(Te,Se) and Fe(Te,S) at the same 
temperature^^ The ratio of J^ oss (H)/j" ±(101) (H) is 
about 5.4 at 1.8 K. Moreover, above 4 K, the critical cur- 
rent densities decrease with the applied field more quickly 
than below 4 K, suggesting that the pinning mechanism 
may change at 4 K. 



IV. CONCLUSION 

In summary, we present the superconducting proper- 
ties of /3-FeSe single crystals. The results indicate that 
the H C 2(T) is isotropic and the spin-paramagnetic ef- 
fect may be the dominant pair-breaking mechanism for 
H||(101) and H_L(101). The calculated GL parameter k 
indicates that /3-FeSe is a typical type-II superconductor 
with the large k. The resistivity exhibits clear Arrhe- 
nius TAFF behavior with a crossover from single vortex 
pinning region to collective creep region for both field di- 
rections. On the other hand, the critical current density 
for field along (101) plane is about five times larger than 
for field normal to (101) plane but still smaller than in 
Fe(Te,Se) and Fe(Te,S) at the same temperature. 
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